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ABSTRACT 

We study the variability of the warm absorber and the gas responsible for the emission lines in the Seyfert 1 galaxy NGC 5548, in 
order to constrain the location and physical properties of these components. Using X-ray spectra taken with the Chandra-hETGS 
in 2002 and 2005, we study variability in the ionic column densities and line intensities. We find a lower Ovn forbidden emission 
line flux in 2005, while the Fe Ka line flux stays constant. The warm absorber is less ionized in 2005, allowing us to constrain its 
location to within 7 pc of the central source. Using both the observed variability and the limit on the FWHM of the O vii f line, we 
have constrained the location of the narrow line region to a distance of 1 pc from the central source. The apparent lack of variability 
of the Fe Ka line flux does not allow for a unique explanation. 
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1. Introduction 

Active Galactic Nuclei (AGN) have been studied for decades 
since their discovery over 60 years ago. Mass loss from the nu- 
cleus has been kn own for a lon g time to exist in the form of 
radio jets (see e.g. Peterson 1997) or in outflows showing broad 

■ absorption lines (.Wevmann et al...l981.) . It is only recently how- 
, ever that X-ray and UV observations have shown the presence of 

outflows in a majority of the moderately luminous Seyfert galax- 
ies (Crenshaw et al. 1999). These outflows are important as they 
probe the inner regions of AGN and during the li fetime of an 

■ AGN they can constitute a significant mass loss ( Blustin et al.l 
[ l2005h . They can provide us with a better understanding of 

the accretion process onto the supermassive black hole and of 
the en richment of the intergalactic medium (Hamann & Ferland 
I Il999l) . The physical structure of these outflows and their con- 
nection to the other AGN components, such as the supermas- 
sive black hole, the accretion disk and the broad and narrow line 
region (BLR, NLR), is not understood. In order to learn more 
about these outflows and the other components of the AGN, vari- 
ability studies are very important. These allow one to see how the 
warm absorber and the emission lines respond to a variable flux 
from the central source. This puts constraints on the locations of 
the gas in which they originate, and provides us with information 
on their physical state. 

A good object for the study of these processes is the Seyfert 
1 galaxy NGC 554 8. It is a relatively nearby AGN (z = 0.01676, 
ICrenshaw & Kraem er 1999) and has a relatively high X-ray 
brightness. Previous UV and X-ray observatio ns have shown the 
prese i ice of a warm absorbe r in NGC 5548 (.Steenbrugge et al.l 
120051; ICrenshaw et al.l 120031) . From the kinematics it was con- 
cluded that the UV and X-ray absorption lines are part of the 
same outflow. By comparing a Chandra-HETGS and LETGS 



observation of NGC 554 8 in 2002 with an earlier one in 1999, 
ISteenbrugge et alj ( 12005 ') conclude that there is no evidence for 
long-term variabihty, except for the O v absorption line. 

We present the results of a followup Chandra-LETGS ob- 
servation of NGC 5548 taken in April 2005, which we compare 
with prev ious Chandra-LETGS and XMM-Newton RGS obser- 
vations (Kaas tra et alll2002bl: ISteenbrugge et aI.l l2003L l20()5). In 
all previously mentioned high-resolution observations, the vari- 
ability in the continuum flux has been at most a factor of 2 and no 
significant variations have been detected in the emission lines. 
During our new observation the source was at a very low flux 
level which lasted at least from March 1 till May 10 , as de- 
termined fro m Swift observatio ns dGoad & Pagell20()6l) . optical 
observations (jBentz et al.ll2007l) . and this LETGS observation. 
During this time the continuum flux was a factor five lower than 
in 2002. This allows us to analyze the source for the first time at 
a very low flux level, and determine the changes in the spectral 
features, such as the warm absorber, naiTow emission lines, and 
the iron Fe Ka line. The Rossi X-ray Timing Explorer (RXTE) 
monitoring data of NGC 5548 from 1996 till 2007 are also pre- 
sented in order to constrain the long-term variability much more 
accurately than only using the spectral observations, which have 
large gaps between them. We discuss the observation and data 
reduction in Section |2] Section [3] is used to describe the spec- 
tral analysis and the variability of the spectral components. We 
discuss our results in Section |4] and present our conclusions in 
Section|5] 



2. Data reduction 
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NGC 5548 was observed for 141 ks with the Chandra LETGS 
(Brinkman et al. 2000) in April 200 5. The data were reduced as 
described bv lKaastra et alj (l2002bl) . In summary, the data were 
reduced with the standard CXC pipeline up until the level 1.5 
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Fig. 1. The rebinned spectrum of NGC 5548 as taken with the 
LETGS instrument in 1999 (upper) 2002 (middle) and 2005 
(lower). 



Table 1. Oxygen column densities for the 2005 observation. The 
logarithms of the column densities are listed in m"^. The veloc- 
ity broadening cr and the outflow velocity u are listed in km s"' . 
The lower table is for the same slab model, bu t with the veloc- 
ity bro adening cr fixed to the vales used in .Steenbrugge et al.l 
(I2005h . 



V 

a 
Ion 



-1040 -530 -160 

127 ±28 62 ± 21 120 + 34 



total 



Oiii 


<20.5 


<20.6 


<21.2 


<21.4 


Oiv 


20.8 ± 0.5 


<20.9 


21.0 ±0.5 


21.2 ±0.5 


Ov 


20.9 ± 0.6 


<21.3 


21.4 ±0.5 


21.5 ±0.5 


Ovi 


22.1 ±0.2 


<20.9 


<21.4 


22.1 ±0.2 


Ovii 


<21.0 


<20.4 


21.1 ±0.5 


<21.6 


VIII 


<21.1 


<21.7 


21.2 ±0.5 


<21.9 


cr 


40 


100 


90 




Oiii 


<20.4 


<20.4 


<21.1 


<21.2 


Oiv 


20.6 ± 0.6 


<21.0 


21.1 ±0.4 


21.3 ±0.3 


Ov 


21.3 ±0.4 


20.5 ± 0.7 


21.4 ±0.3 


21.7 ±0.2 


Ovi 


21.9 ±0.2 


<21.3 


21.1 ±0.4 


22.0 ± 0.2 


Ovii 


21.5 ±0.9 


<20.6 


21.0 ±0.5 


21.6 ±0.5 


VIII 


<21.5 


<21.5 


21.1 ±0.9 


21.3 ±0.6 



event files. For the steps leading to the final event file, level 2, 
we followed an independent procedure, which is described in 
the aforementioned article. The resolution of the LETGS is 0.05 
A and its band spans a wavelength range of 1 - 180 A. Because 
of the low signal to noise ratio and the Galactic absorption to- 
wards NGC 5548, we ignore the data above 80 A and below 1 
A. The spectrum was a nalyzed with the SPEX software pack- 
age (lKaastraetal.lll996fPl . We use Hq = 70 km s' Mpc', Qn, 
= 0.7 and Qa = 0.3. Because of the low count rate we calcu- 
late the errors using C-statistics (see XSPEC online manuaQ, we 
use the form as given by Castor See also Cash 197 9). All errors 
given are lo" errors, co mputed for AC - 1 . The RTXE data were 
reduced as described in lMcHardv et alj (|2004|) . 



3. Spectral Features 

Fig. [T] shows the spectra from the 1999, 2002 and 2005 obser- 
vations. The flux levels are much lower in the most recent data, 
with a 2 - 10 keV continuum flux of 1.47 x lO"'"* W mT^. To 
make sure that the spectral lines are analyzed against an ade- 
quate local continuum, we fit the continuum with a spline model. 
A spline model allows us to model the continuum between two 
boundaries b\ and b2 with a cubic spline. The limits are 1 A and 
80 A in our case. The continuum flux yi is then given on 159 
evenly spaced grid points xi between b\ and ^2 with a spacing 
of 0.5 A between the grid points. The continuum model between 
the grid points is then determined by cubic spline interpolation. 
The y values of the grid points are determined by spectral fitting. 
Using a spline model means we do not assume an a priori shape 
of the continuum, such as for example a powerlaw and black- 
body continuum. We apply the redshift of 0.01676 for NGC 
5548 and correct for Galactic extinction, with the Hi column 
density fixed to 1.65 x 10^"* mT^ (N andr aet al. 1993). We add 
various absorption components and narrow emission features to 
this continuum (see Sect 3.1 - 3.3). The best fit (including all 
the features) has C - 3228 for 2833 degrees of freedom. 



' See also http://www.sron.nI/spe^ 

^ http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/ 



XSappendixCash.html 



3.1. Warm Absorber 

In order t o model the warm a bsorber we use the slab model 
in SPEX JKaastra et al.ll2002ar) . which calculates the transmis- 
sion of a thin, irradiated slab of matter with a set of ad- 
justable ionic column densities. We model the optical depth 
of the absorption lines with Voigt profiles. Our models consist 
of 3 slab components, with the ionic column densities set free 
and w ith the outflow velocities v fixed to the previous LETGS 
study dSteenbrugge et aL 2005h. Simu ltaneous UV observations 
of NGC 5548 (Crensh aw et'ai1l2003h show five different warm 
absorber outflow velocities, but due to lower spectral resolution 
in the X-ray, only three outfl ow velocities could be distinguished 
bv lSteenbrugge et al.l(l2005h . namely -1040, -530 and -160 km 
s-i. 

A change in the optical depth of an absorption line can be 
due to a change in the velocity broadening cr or a change in the 
ionic column density A^ion. Alternatvely a change in the covering 
factor of the absorbing gas can also explain the variations in the 
absorption lines. We assume here a covering factor of unity for 
the absorbing gas (see Sect|4]for discussion). To be certain that 
the changes we observe in the absorption lines are indeed due to 
column density changes and not only due to a possible change in 
cr, we performed two fits, one with cr set free, th e other with cr 
fixed to the values reported in ISteenbrugge et al] (l2005h . Fig. |2] 
shows the 2002 and 2005 total column densities for both models. 
The O VIII column density has decreased, while the O iv - O vi 
column densities have increased. The total column densities are 
the same for both models. In Table[T]we show the 2005 oxygen 
column densities for both models. 



3.2. Narrow emission lines 

The O VII forbidden line at 22. 101 A is clearly visible in the spec- 
trum of the 2005 observation. We also searched for the Ne ix 
forbidden line at 13.699 A, since that line is also detected in 
earlier observations. Table |2] shows the strengths of both lines 
in the new spectrum and in earlier observations. We measure 
a flux of 0.35 + 0.06 ph m"^ s ' for the Ovii forbidden emis- 
sion line and we detect a blueshift of -250 ± 70 km s"' for the 
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Table 2. The unabsorbed flux in ph nr^ s"^ for the O vii and 
Ne IX forbidden emission Unes for observations taken in 2000, 
2001 (RGS, extracted from the pubHc archive, errors include 
uncertainty due to continuum level) and in 2005. For compar- 
ison we also give values from earlier publications for observa- 
tions taken in 1999 (LETGS, Kaasti-a et al. 200 23), 2000 (MEG, 
Kaastra et al.1 l2002 b) and 2002 (LETGS, iSteenbrup^e et al.l 



2005h . All errors are calculated with 68 % confidence. 



Year (month) Instrument 



Ovn 



Neix 



1999 (December) 


LETGS 


0.81 ±0.13 


0.25 ± 0.07 


2000 (Februai-y) 


MEG 


0.82 + 0.18 


0.09 ± 0.03 


2000 (December) 


RGS 


0.65 ±0.18 


0.14 ±0.07 


2001 (July) 


RGS 


0.70 ±0.16 


0.22 ± 0.07 


2002 (January) 


LETGS 


0.88 ± 0.08 


0.14 ±0.04 


2005 (April) 


LETGS 


0.35 ± 0.06 


<0.04 







/, 
/\ 



Fig. 2. Total oxygen column densities as measured with the 

LETGS instrument in 2002 (rectangles) and 2005 (diamonds). le ib 2 22 24 

^ , jj p. 1.111 Wavelength (A) 

Icr errors are shown as well, except tor ions which only have 

a measured upper limit on the total column density (arrows). Fig. 3. The Fe Ka emission line fit for the 2005 observation. 
The upper panel shows the total column densities with the ve- 

lo city broadening cr s et free , the lower for cr fixed to the values Table 3. Fe Ka line parameters from the HEG spectrum (as re- 
oflSteenbrugge et al.| (|2005|). ported by Yaaoobetal. 2001), XMM -Newton EPIC spectrum 

(iPoundset al.1 120031) . HEG spectrum JSteenbnigge et al.ll2005l) 
and present LETGS spectrum (2005). The FWHM is given in 
O vn f line. For the Ne ix forbidden emission line we measure km s' , the flux in ph m^ s' . All errors are calculated with 68 
an upper limit of 0.04 ph m"^ s"', if we assume that the Neix % confidence, 
line has the same blueshift as the Ovn f line. Earlier observa- 
tions as noted in Section [3] have reported blueshifts of -150 + 
70 km S-' (Steenbruggeetal. 2005. HEG) and -70 + 100 km 
s~' dKaastra et al.ll2002bl) . The LETGS has a wavelength scale 
uncertainty of 0.01 A, or 140 km s~' at the O vn f emission line. 
Taking this and the error on the outflow velocity into account, 
the blueshift we detect is negligible, but it is however c onsistent 

with the HEG measurement ( Steenbrugge et a lj|200l). The de- observations. There is no evidence for a relativistically broad- 
crease in flux with respect to previous observations is consistent ^^^^ p^ ^^ emission line profile in our spectrum. We measure a 
for both lines, although we only detect an upper limit for the pjj^^ ^^ 73Q0 + 5 iqo km s' , which is also in agreement with 
Ne IX f hne. previous values. 



HEG(2000) EPIC(2001) HEG(2002) LETGS(2005) 


E (keV) 


6.40±0.02 6.39±0.01 6.39±0.01 6.39±0.03 


EW (eV) 


130±30 60±9 47±11 420±140 


FWHM 


4500±1600 6500 ± 1300 1700 ± 1100 7300±5100 


Flux 


0.36±0.10 0.38±0.06 0.24±0.06 0.55±0.18 



3.3. The narrow Fe Ka line 

The low continuum flux facilitates the detection of a clear Fe 
Ka emission line at 6.39 + 0.03 keV in the spectrum of our 
LETGS data (Fig.[3]|. Using a Gaussian fit, we measure a flux of 
0.55 ± 0.1 8 ph nr^ s"', which agrees with previous valu es from 
lYaaoob e t al. (2001), Pounds et al. (2003) and Steenbrug ge et alj 
( l2005h as shown in Table [3] The equivalent width of the line 
is 420 + 140 eV, much larger than measured in those previous 



4. Discussion 

4.1. The location of the warm absorber 

Fig.|2]shows that the O viii column density has decreased, while 
the O IV - O VI column densities have increased between 2002 
and 2005, both for cr fixed and cr set free. Some errors on the 
oxygen column densities in the 2005 spectrum are large, due 
to a combination of low statistics and saturation of some of the 
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22 22.5 

Restframe wavelength (A) 

Fig. 4. Comparison of the 1999 (diamonds), 2002 (solid line) 
and 2005 (rectangles) spectra, near the forbidden Ovii emis- 
sion line at 22.101 A. The deep Ov absorption line is visible 
at 22.354 A. The 1999 and 2005 spectra have been scaled to the 
2002 flux. 



absorption lines. As shown in Fig. |4] there are additional indi- 
cations that the column density of the lower oxygen ions has 
increased from 2002 to 2005. The O v absorption line at 22.354 
A is indeed much stronger and broader than in either 1999 or 
2002. The equivalent width of the absorption line is 140 ± 30 
mA, compared to 44 + 6 mA in 2002 and 22 + 9 mA in 1999. 
Allowing the velocity width cr to vary, the column densities of 
the individual velocity components do change, but the total col- 
umn densities do not change significantly (see Table [T]i. This is 
a further indication that the observed changes in the absorption 
lines are indeed due to variations in the column density and not 
only due to a possible change in cr. 

In principle this still leaves open the possibility of a vary- 
ing covering factor of the x-ray absorbi ng gas. As is shown i n 
UV spectral observations of NGC 5548 (ICrenshaw et al.ll2003h . 
the covering factor for the various components is not unity. In 
the UV one can resolve the diff'erent velocity components of the 
absorber, while this is very hard, if not impossible to achieve 
with the resolution of the current X-ray instruments. The size 
of the X-ray emitting region (10 Rg is estimated to be much 
smaller than the size of the UV emitting region. Also the BLR 
region is larger for the optical / UV than for the X-rays (see e.g. 
ICostantini et al.ll2007h . This, along with the fact that the column 
densities for the highly ionized oxygen species have decreased, 
while column densities of the lowly ionized species have in- 
creased, allows us to be quite confident with the assumption of a 
covering factor of unity for the X-ray absorbing gas. 
The large increase in O vi can be explained by recombination of 
O VIII and O vii. 

The ionization parameter ^ is deter mined by the ratio of the 
ionizing flux and the density of the gas (lTarteretalJll969l) : 



^ nR^' 



(1) 



Here L is the 1 - 1000 Rydberg luminosity, n the density of the 
gas and r the distance from the ionizing source. In fact, one can 
only determine the product of the density and distance squared, 
nR^, from the luminosity and ^. The luminosity L - 9.0 x lO"*^ W 



in the 2002 LETGS observation dSteenbrugge et al.ll2005h . We 
will first consider the changes in the O viii column density to ob- 
tain a lower limit for the density n. The column density of O viii 
in flie 2002 LETGS ob servation is 3.39 x lO^^ (See model B, 
ISteenbrugge et alj|2005l) . We assume that the value for the ion- 
ization parameter is the one where the column density of the ion 
peaks, so log ^ = 1 .65 for O viii, using the XSTAR codqjfor pho- 
toionization modelling. We assume that the changes we detect 
are due to recombination to determine n. The density n scales in- 
versely with the recombination time Tree of the gas and depends 
on the recombination rat e of the specific ion JKrolik & KrissI 
ll995l:lBottorff^etal.ll2000l) : 



T,-ec(^,) = kr(^/)« 



fiX^) a,{Xi^i) 



f(Xi) a,(Xi) 



(2) 



where a^iXi) is the recombination rate from ion X,=i to ion X, 
and f{Xi) is the fraction of element X in ionization state /. If Tree 
is positive, it means that the amount of O viii increases, by re- 
combination of O IX. If Tree is negative, then this means that more 
O VIII is destroyed by recombination to O vii than that which is 
created by recombination of O ix. So the sign of Tree determines, 
whether O viii is destroyed or created and the value of Tree shows 
how long it takes before equilibrium is reached. 

As the change in column densities occurred between Jan 21 
2002 and April 15 2005, the upper limit on Tree is 1 160 days (in 
the rest frame of the source). The /(X,+i) to /(X,) for i - O viii 
ratio can be determined from the ionization parameter ^. For log 
^ = 1 .65, the fraction of O ix to O viii is 0.48. From ^ we derive 
a temperature of 86 000 K. With this value of T, a^iO vii) = 
8.65x lO'"^m^s-' andar(Oviii)= 1.31 xlO^'^m^s-' dNahad 
1 19991) . Using these values we find a lower limit on the density 
of 4.2 X 10^ m"-'. Using this limit for n and using ([TJ, we get an 
upper limit of the distance of R <1 pc. If we instead use O v to 
determine the distance, using log ^ = -0.2 and T - 19 000 K, we 
get an upper limit of R < 19 pc 

Another method of determining the location of the warm ab- 
sorber is to use its column density Nh and its thickness AR, 
namely Nh - n x AR. We again use the O viii column densit y 
(log f = 1.65, Nh = 7.6 x lO^^ m^ ISteenbrugge et al.ll2005h . 
Since the gas has responded to the lower flux within 1 160 days, 
the upper limit on AR is the light travel distance of 1 pc. From 
this we get a lower limit on n of 2.5 x 10^ m"-'. Again using ([TJ, 
we get an upper limit of the distance of R < 9 pc. 

The main uncertainty causing the large upper limits is of 
course the large gap in the spectral data between 2002 and 2005. 

4.2. The variable O vii f narrow emission iine 



The Une flux of the O vn f line has decreased from 0.88 + 0. 12 to 
0.35 + 0.06 ph m-2 s-> between 2002 and 2005 (Table©. This 
immediately gives us an upper limit to the distance, namely the 
light travel distance of 1 pc. A lower limit of the distance to NLR 
gas can be obtained from the velocity di spersion and th e mass of 
the supermassive black hole as given bv lNetzen (Il990l) : 



M = 



3ry^ 

-" '^FWHM 

4G ' 



(3) 



with M the mass of the supermassive black hole, Vfwhm the 
FWHM of the line and r the distance of the line emitting re- 
gion. This assumes an isotropic velocity distribution, which is 



http://heasarc.gsfc.nasa.gov/xstar/xstar.html 
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reasonable, since we can not obtain any detailed line profile in- 
formation for the NLR region. Using M = (6.54 + 0.26) x 10' 
Mq (Bentz et al. 2007), and an upper limit to Vfwhm of 560 km 
s"', which we have obtained by using a Gaussian line model for 
the O VII f line instead of a delta line, we obtain a lower limit of 
1 .2 pc for the location of the NLR. If instead the NLR would 
foUow the same dist ribution in velocity as the BLR in AGN 
dPeterson et al.ll2004l) . then the lower limi t would be 0.9 pc. A 
similar limit of ~ 1 pc was also derived bv JKaastraet al.l (l2003l) . 
Although this lower limit can change by a factor of a few, de- 
pending on the exact velocity distribution and the upper limit 
can be in fact smaller as it is based on sparse temporal sampling, 
both limits used together constrain the location of the NLR to 
about 1 pc. This is the first time that a well-constrained location 
has been determined, due to the combination of variability and 
line width limits. Optical measurements ( Krae mer et al.l [l998l) 
indicate that the optical narrow line emission comes from a re- 
gion within 70 pc, with a high ionization component at 1 pc, in 
agreement with our findings. 

We have checked whether Tree is much smaller than rit, so 
that the total delay is determined only by rit: 



'''delay - ^It + ^rs 



(4) 



with Tdeiay the total delay between changes in the continuum flux 
and changes in the emission line flux, rit the light travel time 
through the emitting region and r^e the recombination timescale. 
Using ([T]i, with log ^ - I (where the O vii fraction in the gas 
peaks) and R - I pc, we get a lower limit to n of 1 xlO'^ m"-'. 
Again from ^ we can determine the temperature for the gas with 
log ^ = 1, which is 33000 K. The fraction O viii / O vii = 0.25 for 

^ = "l- 

Using the O vii and O vi recombination rates from iNahad 

(IT999L a,(Ovu) = 1.65 xlO"'' m^^s"' and ^^(Ovi) = 0.98 
xl0~'' m^s~'), the upper limit on Tjec is 3 days, much less than 
the 1160 days limit from variability arguments. This means the 
total delay time rjeiay is determined only by the light travel time 

Tit. 

4.3. The origin of the Fe Ka line 

Due to the low continuum flux, we are able to detect the Fe Ka 
emission line with the LETGS, even though the LETGS is not 
optimized for the wavelength range where the line is present, 
since its spectral resolution is 50 mA. The interesting result is 
that the line flux is in agreement with earlier observations, while 
one would expect a decreased line flux, if the line flux correlates 
with the continuum flux. Due to the low continuum flux and the 
constant line flux the equivalent width (EW) of the line is large, 
420 + 140 eV. 

Using both HEG measurements of the Fe Ka line (Table [3]), 
the weighted average of the FWHM is 2600 + 900 km s ' . Once 
again using (O, the line is formed at 0.06 + 0.04 pc from the 
central source, assuming again a mass of ~ 6.54 x 10' Mq. At 
this distance one would expect that the iron line would also have 
responded to the lower ionizing flux, since the light travel time 
is much less than 3 years. 

There are at least five possible explanations for the constant 
flux and corresponding larger EW of the emission line. The first 
one which can immediately be ruled out is the light-bending near 
the black hole (Miniutti & Fabian 2004). The derived distance 
from the FWHM of the line (0.06 pc, or 19 000 Rg), is much too 
large for any light-bending to play a role. 



The second reason why the iron line is not affected by the 
lower continuum could be that the line consists of multiple unre- 
solved components originating from a broad range of distances 
from the central source. If we assume that this is the case, then at 
least part of the emission could come from a region further away 
from the central source than 1 pc, for instance the torus structure 

1 r — 1 

at a parsec-scale (Antonucci 1993). However Sug anuma et al.l 
(2006) have performed infrared reverberation mapping of NGC 
5548 and found a time-lag of 47-53 days between the contin- 
uum flux in the V band and the IR emission (dust). This would 
put the dust at ~ 0.04 pc. The dust can of course extent much 
further outward, so Fe Ka emission from reflection on dust fur- 
ther than 1 pc away from the central source cannot be ruled out. 
However, if the bulk of the Fe Ka line would originate from pc 
distances, the line should be much narrower than what is ob- 
served. 

Another explanation is that the continuum experienced a 
peak in flux just before our observation, such that the Fe Ka 
line would have already responded to this increase in continuum 
flux, but the warm absorber and the narrow emission lines did 
not. From the LETGS observation alone this cannot be ruled out, 
but there were als o optical observations (March 1 till April 10, 
iBentz et al.l[2007h of NGC 5548 just before the LETGS obser- 
vation. They detected a very low continuum and also a weaker 
HyS emission line. Since the H/3 line has a time-lag of 6.3 days, 
any peak in continuum flux should have shown up in the optical 
observations. Also the RXTE monitoring does not show a peak 
in the months prior to our observation. 

The fourth possibility is based on the change in continuum 
flux and the effect it has on the gas which reflects the Fe Ka line. 
Matt et al. (1993) have shown that the strength of the Fe Ka line 
depends on the ionization parameter ^ of the reflecting gas. They 
show that if ^ < 100, the line strength increases with decreasing 
^. Since the continuum flux has dropped, it is not unreasonable 
to assume that the ionization state of the gas producing the Fe 
Kff line will be lower So this means that the line flux could stay 
constant when the continuum flux drops. 

The final possibility is that while the continuum flux in the 2 
- 10 keV band has gone down dramatically, the reflection com- 
ponent above 10 keV which produces the Fe Ka line may be 
much less variable than the continuu m. This has been de tected 
recently in several AGN with Suzaku dReeves et al.ll2006r) . The 
reason for the lack of strong variability in the reflection compo- 
nent is still unknown. Due to the errors on the Fe Ka line flux or 
the instrumental resolution of the LETGS, we cannot distinguish 
among these explanations. 

4.4. Constraints from RXTE monitoring 

Figure |5] shows the long time scale light curve of NGC 5548 as 
observed with RXTE from 1996 till 2007. The data have been 
re-binned to a time scale of one month. As can be seen, the vari- 
ability on a time scale of a few months can be larger than that of 
longer time scales, i.e. years. There still are long-term changes 
visible. From 2002 on, the average count rate is 20 - 25 % lower 
than that of the earlier years. In principle with the information 
we have on the continuum flux of NGC 5548 between 2002 and 
2005, we can further refine our estimate of the upper limit on the 
recombination time scale for the warm absorber 

For the upper limit to the variability time scale of the warm 
absorber we currently use an upper limit of 1160 days. Based 
only upon the two LETGS spectra, this timescale is wildly un- 
certain. We note that the RXTE continuum flux was at approx- 
imately the same level in mid-2004 as it was during the 2002 
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Fig. 5. RXTE lightcurve of NGC 5548, re-binned to one month 
intervals, spanning from 1996 till 2007. The count rate shown is 
per PCU detector of the RXTE. All 6 spectral observations (see 
Table 121 are shown as well. There is a gap in the data from the 
beginning of 2003 till 2004, when the source was not monitored 
with RXTE. The errors on the countrates are at the 1 % level. 



LETGS observation. So if the ionization timescale is short, the 
ionization state in mid-2004 might be the same as in 2002. In 
which case a better estimate for the recombination timescale 
might be 1 year, and even that might be far too long. E.g. there 
is a peak in the X-ray flux only a few months before our 2005 
observation, so an upper limit of a few months might be even a 
better estimate. If the recombination timescale is 4 months in- 
stead of 1160 days, the upper limit on the warm absorber dis- 
tance would be 2 pc instead. As this recombination timescale is 
so uncertain, and likely to be a serious over-estimate, we there- 
fore conclude that the warm absorber is probably much further 
in than the 7 pc we derived as an upper limit. Further studies, 
taking time-dependent ionization into account, will be able to 
better model the changes in the warm absorber as it reacts to the 
continuum variations that RXTE has observed. 

Concerning the changes in the emission lines, previous ob- 
servations have shown no observable change in flux for both 
lines (see Tables|2]and[3]), even though the continuum varied sub- 
stantially between those observations. This is a clear indication 
that both emission lines respond to the continuum changes on a 
long time scale, i.e. years not weeks or months. For the Fe Ka 
line this would mean that any correlation between the continuum 
flux and the line flux is more complex than a linear relation be- 
tween the two. Possible correlations between the emission lines 
and continuum flux will be investigated further in a future paper. 

5. Conclusions 

This is the first detailed X-ray analysis of NGC 5548 at an ex- 
tremely low flux level. The warm absorber which has shown lit- 
tle signs of change in earlier observations shows clear signs of 
recombination to a lower ionized state. We have constrained its 
location to within 7 pc of the central source, although this is most 
likely a serious overestimate. The narrow forbidden O vii line has 
decreased in flux, which along with the limit on the line width, 
places its location at ~ 1 pc from the central source. This loca- 
tion also agrees with previous estimates for this source. Only the 
Fe Ka emission line has not responded linearly to the lower con- 



tinuum flux, the line flux is consistent with earlier observations. 
The most likely explanations are a change in the ionization state 
of the reflecting gas or a constant reflection component which 
does not follow the 2-10 keV continuum. 
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